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Abstract 


We derive an expression connecting the coefficients of a polynomial ex- 
panded in the Bernstein basis to the coefficients of an equivalent ex- 
pansion of the polynomial under an affine mapping of the domain. The 
expression may be useful in the calculation of bounds for multi-variate 
polynomials. 


1 Introduction 


The set of Bernstein basis polynomials of degree n can be used to form 
a basis for a vector space of polynomials of degree less than or equal to 
n. These polynomials take the form 


Bi 


n 

v 


x u (l — x) n ~ v , is = 0, . . . ,n 


(1) 


where 



is a binomial coefficent. It is known that the coefficients 


of a given polynomial expressed using the Bernstein basis provide guar- 
anteed bounds on the global minimum and maximum of the polyno- 
mial [1,2]. This result extends from the fact that these basis polynomials 
are non-negative over the unit interval and the values at the end points 
are simply 

Bv,n(fi) — <U and B vn (X) — Tm (2) 


where 5 is the Kronecker delta function. Furthermore, such bounds may 
be improved by calculating and comparing the various expansion, or 
Bernstein, coefficients for the polynomial over subdivisions of the do- 
main interval, i.e. the coefficients that result from an affine mapping of 
a subdivision of the domain back to the original domain. The de Castel- 
jau algorithm [3] is a method for calculating the Bernstein coefficients 
that are generated via a repeated halving of the domain. Munoz and 
Narkawicz [4] have provided a simpler variant of this algorithm, by ex- 
plicitly calculating the analytical expressions that relate the coefficients 
on the original and halved domains. We present here an expression con- 
necting the Bernstein coefficients on the unit domain to those for an 
affine mapping of the domain. This is in a sense a generalization of the 
formulas given by Munoz and Narkawicz [4], The formula here reduces 
to their result in the special case of dividing the domain into halves. We 
also show how to apply the equation to a division of the domain into any 
number of intervals of equal size. 


2 Derivation of the General Formula 

We derive here the relationship between the coefficients of a single vari- 
able polynomial in Bernstein form on the unit interval x £ [0, 1] , B UjTI (x), 


1 



and those of the corresponding Bernstein polynomial generated by the 
mapping x — > ax + j3. This result holds for the multivariate case as can 
be easily shown using Smith’s representation [5]. The derivation makes 
use of two results: 1) the connection between the coefficients of an ex- 
pansion in B l/tn (ax) to those of an equivalent expansion in B Vin (x ) and 
2) the identity 

B u , n { 1 -x) = B n - Vtn {x) (3) 

The three steps shown below create the mapping x — > ax + /?. The first 
result is used in the first and third steps, and the second result is used 
for the second step. 


x'i = a\x 

(4) 

X 2 = 1 — x\ = 1 — a\x 

(5) 

x 3 = a 2 X 2 = a 2 ( 1 - a\x) 

(6) 


where (5 = a 2 7^ 0 and a = —a\f3. The derivation of the relationship 
proceeds similarly to the procedure in Section 2.3 of Ref. [4], and begins 
by determining the appropriate partition of the argument 1 — ax. 
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v 
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a v x v { 1 — x + (1 — a)x) r 
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k = 0 


n — v 


(1 — a) k x k (l — x) n ~ p ~ k {7) 


n—v—k( 
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k=u 


k 

n — v 
k — v 


(1 — a) k ~ y x k {l — x) n ~ k (8) 


n—k 


Using the trinomial revision formula 


k \ n 


n \ n — i 


we have that 


B u ,n(cxx) = a u ^2 
k=v 


k 

v 


k — i 


(1 - a) fc -"R fc , n (x) 


(9) 


( 10 ) 


Equating two equivalent expansions for the same polynomial, one in the 
Bernstein basis with argument x, and one in the Bernstein basis with 
argument ax, 

n n 

^ ) ('v.nBu.n. (ax) = 'y ) Ck, n Bk,n{x ) (1 1) 

u=0 k = 0 

gives the following identity involving the Bernstein coefficients c“ n and 

k / \ 

*,n = E <» „ “"(1 - = E ( 12 ) 

u= 0 V / v=0 
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where the ainc^ n is not an exponent but a notational superscript. Using 
this result, the identity in Eq. (3), and representing now the Bernstein 
coefficients of an expansion in the arguments x, x' 2 , and X 3 as bk. n , b\ n , 
and b\ n respectively, the coefficients corresponding to Eqs. (4), (5), and 
( 6 ) have the following relationships: 


n n 


E b l,n B vA a 2 x 2 ) = E b l,n B k,n( x 2 ) 

v=0 k= 0 

(13) 

k 

b l,n = E b l,n B u,k( a 2) 

v= 0 

(14) 

n n 

y ( b u ) n B y,n{ x 2 ) y ] b v,n B v,n( 1 — x l) 

u=0 u=0 

(15) 

= b u,n B n-i',n{x l) 

u=0 

(16) 

Tl 

y ) b n—v,n B v,n{x l) 
u=0 

(17) 

— y ) b k J n B k,n{.X ) 

k= 0 

(18) 


The coefficients in the last expression can thus be written as 

k 

b k,n = b n — Utn B u,k( a l) 

v=0 
k n—v 

= E E (19) 

v=0 v '= 0 P 


This equation obviously holds only for (3^0. When /3 = 0, we can 
simply apply Eq. (12) instead. 


3 Specific Cases 

It is straightforward to compute the coefficients for the case where the 
unit domain is divided into halves. This division is equivalent to the 
mapping x — > ax + (3, where a = 1 / 2 , f3 = 0 , for the left half, and 
a = —1/2, (3 = 1 for the right half. In the former case, we use Eq.(12) 



and in the latter, Eq.(19) yields 
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( 21 ) 


k n—u ^ 

bk,n = FI F F , n Bv,k{-^)B v i >n - u (l) 

u =0 u '= 0 


1 

FI 9 ) 

u=0 

( 22 ) 

Iff k ) b * 

2 fe h\ v ) n ~ u,n 

(23) 


which are the same formulas given in Eq.(14) of Ref. [4]. 

These equations can be used to calculate the Bernstein coefficients 
corresponding to domains that are created by repeated halving of the 
unit domain. Alternatively, we may from the outset define a linear divi- 
sion of the domain into M equal size windows by simply setting 

a = —1/M, j3 = m/M, m = 1, . . . , M (24) 

Such a scheme allows one to repeatedly divide the unit domain into 
smaller pieces faster and perhaps in a non-even fashion if one so chooses. 

4 Conclusions 

We have derived an expression relating the Bernstein coefficients of a 
polynomial on the unit domain x € [ 0 , 1 ] to those for the mapped do- 
main x — > ax + (3. The expression provides a generalization of the 
method implemented by Munoz and Narkawicz [4]. The full implica- 
tions of this work are still being determined; for example, one may study 
the aforementioned linear division of the domain in the limit of infinitely 
small divisions, M — > oo. 
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